o 

(N 

C 

c3 



in: 
> 

OO 
l> 
(N 

in 

(N 

On 
O 



How the recent BABAR data for P — > 77* affect the Standard 
Model predictions for the rare decays P — > 

A.E. Dorokhov a * 

a Joint Institute for Nuclear Research, Bogoliubov Laboratory of Theoretical Physics, 
141980 Dubna, Moscow region, Russian Federation; 

Institute for Theoretical Problems of Microphysics, Moscow State University, 
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Measuring the lepton anomalous magnetic moments {g — 2) and the rare decays of light 
pseudoscalar mesons into lepton pairs P —> l + l~ , serve as important tests of the Standard Model. 
To reduce the theoretical uncertainty in the standard model predictions, the data on the charge 
and transition form factors of the light pseudoscalar mesons play a significant role. Recently, 
new data on the behavior of the transition form factors P —¥ 77* at large momentum transfer 
were supplied by the BABAR collaboration. There are several problems with the theoretical 
interpretation of these data: 1) An unexpectedly slow decrease of the pion transition form factor 
at high momenta, 2) the qualitative difference in the behavior of the pion form factor and the 
77 and rj form factors at high momenta, 3) the inconsistency of the measured ratio of the r\ and 
rf form factors with the predicted one. We comment on the influence of the new BABAR data 
on the rare decay branchings. 
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Modern cosmology tells us that 95% 
of the matter in the universe is not de- 
scribed in terms of the Standard Model 
(SM) matter. New excitement appeared 
after the satellite experiments Fermi LAT, 
PAMELA, ATIC, HESS and WMAP that 
provided data which give an indication 
interpreted as Dark Matter and/or Pul- 
sar signals. Thus the search for hints 
of Physics beyond SM is the fundamen- 
tal problem of modern elementary particle 
physics. There are two strategies to solve 
this problem. 



*The author acknowledges partial support from 
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Firstly, in high energy experiments an 
enormous amount of energy is applied in 
order to excite the heavy degrees of free- 
dom expected to be associated with the 
new physics. At the moment, there is not 
any evidence on a deviation between the 
SM predictions and the high energy data, 
and we are urgently waiting for the results 
of the physical program at LHC. 

Another strategy is to carry out low en- 
ergy experiments, where the lack of en- 
ergy is compensated by a huge statistics 
producing rare processes related to the 
new interactions. The low energy exper- 
iments are not only supplement to the 
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high energy program, allowing to get com- 
bined restrictions on the parameters of hy- 
pothetical interactions, but they are also 
unique instrument for the discovery of the 
physics beyond the SM, containing low 
mass particles. At the moment, there 
are some problems with the matching of 
the experimental data with the predic- 
tions of the SM. The most famous one is 
the descrepancy by three standard devia- 
tions between experiment [1] and SM the- 
ory [2j for the muon anomalous magnetic 
moment (AMM). Another example is sim- 
ilarly large deviation between the recent 
precise experimental result on the rare 7r° 
decay into e~e + pair [3] and the SM pre- 
diction pi5H6lr7] . 

At the early stage of the study of 
the lepton anomalous magnetic moments, 
a i — (9 1 ~~ 2)/2, entering the structure of 
the vector vertex 

ic 

T M = e7 M + at- — a^q,,, (1) 
zmi 

played a fundamental role in the founda- 
tions of quantum mechanics and, in par- 
ticular, of quantum field theory [8|. At 
present, the accent of the study is shifted 
to the test of the SM and the search of 
the physics beyond it. Within the SM, 
the dominant contribution to the lepton 
AMM is due to the QED, supplemented 
by small, but visible corrections from the 
strong and weak interactions. 

The electron AMM is measured with 
one of the best accuracies obtained for 
physical observables [9] 

af p = 1159652180. 73(0.28)-l(T 12 [0.28ppt]. 



In the SM, it is given by 

„SM _ QED 1 hadron , weak 
a e ~ a e ' e a e i 

«. QED = E^(f )" + ..., (3) 

where the first three coefficients are known 
analytically and the two others in some 
approximations. If the fine structure con- 
stant a would be known from other inde- 
pendent sources, the measurement of the 
electron AMM would be a stringent test 
of QED. However, the theoretical error is 
dominated by the uncertainty in the in- 
put value for the QED coupling constant 
a, and the problem is reversed to that of 
obtaining a best estimate of the QED cou- 
pling constant j9] 

a' 1 = 137.035999084(51) [0.37ppb]. (4) 

The great feature of the study of the 
leptonic AMM is, that they are very sen- 
sitive to the manifestation of new physics. 
Any interaction with characteristic scale 
A contributes to the leptonic AMM 
(mi/ A) 2 . Therefore, the heavier the lep- 
ton, the more visible the interaction. In 
this way, the contribution of an interac- 
tion to the muon AMM is bigger than 
that to the electron AMM by a factor 
(m^/rrie) 2 ~ 10 4 . Even bigger would be 
the effect for the r-lepton AMM. How- 
ever, the r-lepton is highly unstable, and 
the measurement of its AMM is problem- 
atic (see [11] for discussions). At moment, 
there are only very rough experimental 

2 An alternative mechanism for a contribution by 
new physics was proposed in [TU]. It occurs at 
the tree level and exhibits a linear rather than 
quadratic dependence on mi. 
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limitations on the r-lepton AMM set by 
the L3 [12], OPAL [KJ] and DELPHI col- 
laborations [mu], 

- 0.052 <a T < 0.058, L3, 
-0.068 <a T < 0.065, OPAL, 
-0.052 <a T < 0.013, DELPHI (5) 

from Z — > TT'-f and e + e~ — > e + e~r + r~ 
processes, while the SM prediction is [11 J 



1.17721(5) ■ 10" 



(6) 



The theoretical studies of the muon 
AMM g — 2 (see for review [T6"lll7|ll8|n9l 
[2D] ) , the rare decays of light pseudoscalar 
mesons into lepton pairs [4"f5f6"fT] and the 
comparison with the experimental results, 
offer an important low-energy tests of the 
SM. The discrepancy between the present 
SM prediction of the muon AMM and its 
experimental determination pQ is (24.6 ± 
8.0) • 10~ 10 (3.1a) [20]. The situation with 
the rare decays of the light pseudoscalar 
mesons into lepton pairs became more 
intriguing after the recent KTeV E799- 
II experiment at FermiLab [3] in which 
the pion decay into an electron-positron 
pair was measured with high accuracy us- 
ing the K L — > 3n process as a source 
of tagged neutral pions (R(P — > l + l~) = 
T(P^l+l-)/T tot ) 

R KTeV ^ e+e -^ = ( 7Ag ± .38)-10- 8 .(7) 
The standard model prediction gives [417] 

^Theor ^0 ^ g+ ^ = (g 2 ± .1).1(T 8 , (8) 

which is 3.1cr below the KTeV result (j7J). 

The main limitation on realistic predic- 
tions for these processes originates from 



the large distance contributions of the 
strong sector of the SM, where perturba- 
tive QCD does not work. In order to di- 
minish the theoretical uncertainties, the 
use of the experimental data on the pion 
charge and transition form factors are of 
crucial importance. The first one, mea- 
sured in e + e~ — > n + n~ (7) by CMD-2 
, SND [22j, KLOE [23], and BABAR 
provides an estimate for the hadron 
vacuum polarization contribution to muon 
o — 2, with accuracy better than 1%. The 
second one, measured in e + e~ — > e + e~P 
for spacelike photons by CELLO [25J, 
CLEO [26J, and BABAR [27] collabora- 
tions and in e + e~ — > P7 for timelike pho- 
tons by the BABAR [28] collaboration, is 
essential to reduce the theoretical uncer- 
tainties in the estimates of the contribu- 
tions of the hadronic light-by-light process 
to the muon o — 2 and in the estimates 
of the decay widths of P — > . The 
BABAR data [28.27J on the large momen- 
tum behavior of the form factors cause the 
following problems for their theoretical in- 
terpretation: 1) An unexpectedly slow de- 
crease of the pion transition form factor 
at high momenta [27J, 2) the qualitative 
difference in the behavior of the pion and 
r), r)' form factors at high momenta [29], 3) 
inconsistency of the measured ratio of the 
77, r)' form factors with the predicted one 
[28]. 

In Figs. 1-3 the data for the 7r°, t] 
and 7/ transition form factors from the 
CELLO, CLEO, and BABAR collabora- 
tions are presented. In Figs. 2 and 3, the 
BABAR point, measured in the timelike 
region [251 . is drawn at Q 2 = 112 GeV 2 , 
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assuming that the spacelike and timelike 
asymptotics of the form factor are equal. 
It is seen from the Figs. 2 and 3, that 
the spacelike and timelike points are con- 
jugated. We hope that new data from the 
BABAR and BELLE collaborations con- 
firm this assumption. 




Figure 1. The transition form factor 7*7 -> ir°. 
The data are from the CELLO [25], CLEO [26] 
and BABAR [27] Collaborations. The dotted line 
is massless QCD asymptotic limit. (The notation 
for curves is explained in the text.) 



At zero momentum transfer, the transi- 
tion form factor is fixed by the two-photon 
decay width 



F 2 



-P77 



.(0,0) 



1 64vrr(P ->■ 77) 



(9) 



where a is the QED coupling constant, 
Mp is the resonance mass and Y{P — > 77) 



> 

CD 

a 



□ CELLO 

CLEO 
♦ BABAR 




Q (Ge\r) 

Figure 2. The transition form factor 7*7 -4 r\. 
The data are from the CELLO [35], and CLEO 
[26 Collaborations. The CLEO results obtained 
in different rj decay modes are averaged. The 
BABAR point, measured in the timelike region 
7* -> 777 [28], is drawn at Q 2 = 112 GeV 2 , as- 
suming that the spacelike and timelike asymp- 
totics of the form factor are equal. (The notation 
for curves is explained in the text.) 



is the two-photon partial width of the me- 
son P. The axial anomaly predicts 



Pit 



(0,0) 



(10) 



where fp is the meson decay constant. 
Under assumption of factorization, per- 
turbative QCD predicts the asymptotic 
behavior of the Fp 77 *(Q 2 ,0) transition 
form factors as Q 2 — > 00 [30] 



F Pir (Q 2 ^oo,0) 



2/f 

Q 



2 ■ 



(ii) 



The perturbative QCD corrections to this 
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CLEO 
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Figure 3. The transition form factor 7*7 -> 77'. 
The data are from the CELLO [25], and CLEO 
[25] [25] Collaborations. The CLEO results ob- 
tained in different 77' decay modes are averaged. 
The BABAR point, measured in the timelike re- 
gion 7* — > 7/7 [25J, is drawn at Q 2 = 112 GeV 2 , 
assuming that the spacelike and timelike asymp- 
totics of the form factor are equal. The dashed 
line is the perturbative QCD asymptotic limit. 
(The notation for curves is explained in the text.) 



expression at large momentum transfer 
are extremely small [31 32 33f34j . 

To describe the soft nonperturbative re- 
gion of Q 2 , a simple interpolation between 
Q 2 — > and Q 2 — > 00 limits has been pro- 
posed by Brodsky and Lepage (BL) 



F^(Q 2 ,0) 



7T77* 



4vr 2 /p 1 



(Ap L ) 2 



^f 2 



pi 



c?7(A? L ) 2 ' 

(12) 



where the values of fp are estimated from 
P and (flO]) [26]: fa = 92.3 MeV, = 
97.5 MeV, /„, = 74.4 MeV. 



The CLEO (and CELLO) collaboration 
parameterized their data by a formula 
similar to (JT21 . but with the pole mass 
being a free fitting parameter [26J , 

F^°(Q 2 ,0) = 47c 2 fpl + Q2 /A 2 p i ( 13 ) 

where = 776 ± 22 MeV, A,, = 774 ± 29 
MeV, and A v > = 859 ± 28 MeV. 

In Figs. 1-3 the asymptotics (fTT!) are 
shown by dotted lines, the BL interpo- 
lations (fl~2l) are given by dot-dot-dashed 
lines, and the CLEO parametrizations ex- 
trapolated to higher momentum transfer 
(fT3"]) are shown by dot-dashed lines. We 
see that the QCD inspired expression f[T2|) 
works well only for the 7] meson form fac- 
tor with A^ L = ^/8ir 2 f 2 ns 866 MeV 
(Fig. 2), whereas the CLEO parametriza- 
tion (fT3|) underestimates the large Q 2 be- 
havior. On the other hand, the CLEO 
parametrization describes the 7/ meson 
form factor (Fig. 3) well, but the BL ex- 
pression strongly underestimates the large 
Q 2 behavior. We still have a good descrip- 
tion of the rj' meson form factor by BL 
formula if one takes = 125 MeV (short 
dashed line in Fig. 3), but then the nor- 
malization is incorrect. 

For the rj and 7/ mesons the 
parametrizations ( TT2|) and (fT3|) correctly 
reflect the experimental data at large Q 2 
on the qualitative level. This is not the 
case for the pion form factor showing the 
growth at large Q 2 , which is unexpected 
from the QCD factorization approach [34J 
(Fig. 1). However, this growth is easy to 
explain [32] in the context of the quark 
model [36]. Within this model, the pion 
form factor is given by the quark-loop 
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(triangle) diagram with momentum inde- 
pendent quark mass serving as an infrared 
regulator [21]. The form factor has dou- 
ble logarithmic asymptotics at large Q 2 : 
ln 2 (Q 2 /M 2 ) and is given by [36] 



•{2 arcsin 2 '' 



ml 



m Tl 



2M C 



ml + Q 2 2 arcsin 2 (^) 

-iln 2 ^±i}. (14) 
2 Bo-V K ' 



where 



4M 2 

1 + -yS-. The solid line in 



Fig. 1 is the pion transition form factor 
calculated from Eq. (fT4j) with the param- 
eter Mq = 135 MeV. The advantage of 
this model is, that it defines the form fac- 
tor for arbitrary virtuality of photons and 
has correct normalization at zero virtuali- 
ties. One of disadvantages is, that the cor- 
responding integral for the decay constant 
f w is divergent and should be regularized. 
Another feature of the model is, that the 
vertices and propagators are local. In par- 
ticular, the pion-to- quarks vertex is lo- 
cal, just like in the Nambu-Jona-Lasinio 
model. It is known, that in this case the 
pion distribution amplitude and distribu- 
tion function are constants [38,39 40 41j. 

The flat (almost constant) pion dis- 
tribution amplitude became popular in 
the context of the explanation of the 
BABAR data within different factoriza- 
tion schemes [4"2l43|44|4"5] . For example, 
in the model [42] the pion transition form 
factor is 

F W ~(~(*(Q ,0) = — 

f 1 dx 
Jo X 




and has logarithmically enhanced asymp- 
totic behavior ~ log (1 + Q 2 /a). In the 
kinematical range of Fig. 1, the model 
(II 5 j) practically coincides with the model 
(HU), if the parameter a = 0.48 GeV 2 . 
Note, that these logarithmically enhanced 
models are not able to describe the 77, 7/ 
form factors. 

However, these approaches still contain 
some unanswered questions. First of all, 
the relation of the parameter a in (|T5|) 
(or m in [43]) to the fundamental QCD 
parameters is unclear. Secondly, the ori- 
gin of the flat distribution amplitude is 
not well justified. Most of the QCD sum 
rule and instanton model calculations lead 
to endpoint suppressed amplitudes (see, 
e.g. [HUT]). This is a simple conse- 
quence of the nonlocal structure of the 
QCD vacuum. The pion-to-quarks ver- 
tex has its own hadronic form factor with 
the characteristic scale of the vacuum non- 
locality (the instanton size). It leads to 
suppression, if realistic values of param- 
eters are used. Only under the assump- 
tion that vacuum nonlocalities disappear, 
a flat pion distribution amplitude is ob- 
tained. Note, that the flatness is a natural 
property of the photon distribution ampli- 
tude, because the photon has no intrinsic 
hadronic form factor [45] . 

The possible origin of the difference of 
the asymptotic behavior of the pion and 
77, 7/ meson form factors is the flavor com- 
position of these mesons |29j. The pion 
consists of almost massless u, d quarks, 
while the 77, 77' mesons include also an s 
quark. The s quark with mass m s of order 
Aqcd m &y be considered as a heavy one. 
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Table 1 



Values of the branchings R (P — > l + l ) — T (P — >■ l + l ) /T tot obtained in our approach 
and compared with the available experimental results. 



R 


Unitary 
bound 


CLEO+BABAR 
bound 


CLEO±BABAR 
±OPE 


With mass 
corrections 


Experiment 


R(n° ->- e+e") x 10 8 


> 4.69 


> 5.85 ±0.03 


6.23 ±0.12 


6.26 


7.49 ± 0.38 [3] 


R(r]^ fj+fT) x 10 6 


> 4.36 


< 6.60 ±0.12 


5.35 ±0.27 


4.76 


5.8 ±0.8 [54 ; 55J 


R(r] e+e") x 10 9 


> 1.78 


> 4.27 ±0.02 


4.53 ±0.09 


5.19 


< 2.7 ■ 10 4 [56J 


R(tf -> x 10 7 


> 1.35 


< 1.44 ±0.01 


1.364 ±0.010 


1.24 




R(i ->• e+e-) x 10 10 


> 0.36 


> 1.121 ±0.004 


1.182 ±0.014 


1.83 





Recently, a similar behavior to that pre- 
dicted by fflUj) was found for the 77* — > r\ c 
transition form factor, measured by the 
BABAR collaboration for the range Q 2 = 
2 — 50 GeV 2 [53]. The corresponding fitted 
mass parameter is A Vc = 2.92 ± 16 GeV. 

Let us consider the possible influence of 
the BABAR data on the rare leptonic de- 
cays of the light pseudoscalar mesons. Re- 
mind, that the imaginary part of the am- 
plitude of these decays is fixed unambigu- 
ously by the two-photon decay widths of 
these mesons. Neglecting the real part of 
the amplitude, the so called unitary bound 
is obtamecjE In particular, one has 



R[7r°^e + e-) (16) 
> # unitary (tt° ->• e+e") = 4.69 ■ 10~ 8 , 

which is 7.3cr below the branching 
R KTeV ^0 e+e -) = ( ? 49 ± 0>38) . 10 - 8 

obtained by the KTeV collaboration [3]. 

The structure of the real part of the 
amplitude was considered in detail in 
[4"fof7] . It consists of the logarithmic 



3 In general the unitary bound related to the two- 
photon intermediate state is violated in the case 
of rf meson where new thresholds appear. 



(model independent) terms ln(m 2 /Mp) 
and ln(mf/A 2 ), a constant term related to 
the inverse moment of the pion transition 
form factor in symmetric kinematics 








-§ ±Xp(m), (17) 



with 



xp^) = — a -- 2 



dt- 



(M) 



+ 



dt- 



7T7*7* 



(t,t) 



t 



and of small mass corrections, (m 2 e /Mp) n , 
(m 2 /A 2 ) n , (M 2 /A 2 ) n , where A « M p is 
characteristic parameter of the form fac- 
tor. As was explained in [I], the data 
on the transition form factors _Fp 77 » (Q 2 , 0) 
provide improved lower bounds for elec- 
tronic decay modes, because in this case 
the logarithmic terms in the amplitude 
dominate over the constant term (see Ta- 
ble 1). On the contrary, for the muonic de- 
cay modes of the 7] and rj' mesons, the log- 
arithmic and constant terms are compara- 
ble, and by using the data on F PjJ * (Q 2 , 0) 
one gets upper bounds for the branchings 
(Table 1). The analysis of the CELLO and 
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CLEO data on the pion transition form 
factor leads to an improved bound 



7r — > e e 



> R CLEO (tt e + e 
= (5.85 ±0.03) ■ 10"' 



which is 4.3a below the KTeV result. 

Let us check how sensitive these im- 
proved bounds are to the recent BABAR 
data at large momentum transfer. As seen 
from Figs. 1-3, there are two main changes 
after the appearance of the BABAR data, 
compared with the CLEO parametriza- 
tion. Firstly, the tail of the pion form 
factor has unexpected asymptotics, and, 
secondly, « 866 MeV is closer to the 
corresponding rf parameter. In order to 
estimate the influence on the pion decay, 
we choose the parametrization suggested 



in 



Fwyj* , 0) 



3.1/f .\/f, 



hQ 4 + h 5 Q 2 + h 7 
(Q 2 + M2)(g2 + M 2 / 



(19) 



where M v = 769 MeV, M v > = 1465 
MeV and the parameter hj = 14.153 is 
fixed by the anomaly (j9l) . The best fit 
to the CELLO and CLEO data is given 
by h s = 6.93 and h x = [49] . The best 
fit, including the BABAR [27] data corre- 
sponds to [50] h$ = 6.51 and switching on 
the small coefficient hi = 0.17, responsible 
for the asymptotics observed by BABAR 
(dashed line in Fig. 1)0. Comparing the 
constant f[T7|) . calculated with one ([TBI or 



4 However, the parametrization (TTO)) can not be 
considered as physical one and used for an extrap- 



another (TT91) parameterization of the data, 
with the integral in (|T7|) . taken in the re- 
gion from to 40 GeV 2 , we find only small 
discrepancy and thus no changes for pre- 
dictions of the 7r~ — >■ e + e~ decay (Table 
1). A similar conclusion was found for 
the influence of the BABAR data [27] on 
the hadronic light-by-light scattering con- 
tribution to the muon anomalous moment 



The new scale A,, for the rj meson form 
factor, obtained from the inclusion of the 
BABAR data [28], slightly changes the 
numbers in the second and third lines of 
the Table. The fifth column of Table con- 
tains the predictions when the mass cor- 
rections to the amplitude are taken into 
account [7J. These corrections have some 
influence on the predictions for the rj and 
rj mesons decays. Thus, it is clear, that 
more precise data at low energy would 
lead to stronger restrictions on the lep- 
tonic decays widths of the light pseu- 
doscalar mesons. 

Further independent experiments for 
tt -> e+e" at WASAatCOSY [57J and for 
rj(r}') -> l+l- KLOE [58] and BES III [59] 
and other facilities will be crucial for res- 
olution of the problem with the rare lep- 
tonic decays of light pseudoscalr mesons. 
It is also important to confirm the theo- 
retical basis for a maximally model inde- 
pendent prediction of the branchings (see 
Table) by getting more precise data on the 
pion transition form factor in asymmet- 
ric as well as in symmetric kinematics in 



olation to higher Q 2 , because it contradicts the 
Terazawa-West inequality [51152] F^ 77 .(<5 2 ,0) < 
l/Q following from unitarity. 
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wider region of momentum transfer. Such 
data are expected soon from the BABAR, 
BELLE (at large momentum transfer) and 
KEDR (at small momentum transfer) col- 
laborations. 

There are quite few attempts in the lit- 
erature, to explain the excess of the ex- 
perimental data on the n° — >■ e + e~ decay 
over the SM prediction, manifesta- 
tion of physics beyond the SM. In Ref. 
[60], it was shown that this excess could 
be explained within the currently popular 
model of light dark matter involving a low 
mass (~ 10 MeV) vector bosons U^, which 
presumably couple to the axial- vector cur- 
rents of quarks and leptons. Another pos- 
sibility was proposed in Ref. [6"Tf62] . in- 
terpreting the same experimental effect as 
the contribution of the light CP-odd Higgs 
boson appearing in the next-to-minimal 
supersymmetric SM. 

The author thanks S.B. Gerasimov, 
V.P. Druzhinin, N.I. Kochelev, E.A. Ku- 
raev, S.V. Mikhailov, A. A. Pivovarov, 
A. V. Radyushkin for discussions on the in- 
terpretation of the high momentum trans- 
fer data for the pseudoscalar meson tran- 
sition form factors. 
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